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ABSTRACT 



f 



We examine the cosmic evolution of a stellar initial mass function (IMF) in galaxies that 
varies with the Jeans mass in the interstellar medium, paying particular attention to the K- 
band stellar mass to light ratio (M/Lk) of present-epoch massive galaxies. We calculate the 
typical Jeans mass using high-resolution hydrodynamic simulations coupled with a fully ra- 
diative model for the ISM, which yields a parameterisation of the IMF characteristic mass 
as a function of galaxy star formation rate (SFR). We then calculate the star formation his- 
tories of galaxies utilising an equilibrium galaxy growth model coupled with constraints on 
the star formation histories set by abundance matching models. Our main result is that at 
early times, energetic coupling between dust and gas drive warm conditions in the ISM, and 
hence bottom-light/top-heavy IMFs associated with large ISM Jeans masses for massive star- 
forming galaxies. At late times, lower cosmic ray fluxes allow for cooler ISM temperatures 
in massive galaxies, and hence bottom-heavy IMFs. Because the massive stars (A/ > IMq) 
formed during the top-heavy phases at early times have all disappeared by today, the resultant 
M/Lk{c') in massive galaxies at the present epoch is increased relative to the non-varying 
IMF case. A key result is that a given galaxy may go through both top-heavy and bottom- 
heavy phases during its lifetime. Quantitatively, the variations in M/Lk{(^) with galaxy mass 
are slightly smaller than what is observed. Nonetheless, these results are encouraging because 
they can, at least qualitatively, reconcile a bottom-light IMF that would help explain a num- 
ber of observations of massive high-z galaxies with the bottom-heavy IMF inferred for the 
descendants of those galaxies today. 

Key words: stars:luminosity function, mass function - stars: formation - galaxies; formation 
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1 INTRODUCTION 

The stellar initial mass function (IMF) is fundamental for under- 
standing a wide range of problems in astrophysics. The IMF pro- 
vides a mechanism for parameterising the distribution of masses 
of stars formed in a given generation. Beyond being a cru- 
cial ingredient for any successful theory of star formation (e.g. 
iMcKee & Ostrikej|2007l) . the IMF determines such galactic prop- 
erties as the metal enrichment history, energy distribution of pho- 
tons injected into the interstellar medium (ISM), and conversion 
of luminosities into physical quantities such as star formation rates 
(SFRs) and stellar masses (Af*). 

Despite the obvious importance of the IMF, no consensus ex- 
ists that describes its origin and possible variations. In the Milky 
Way, direct star counts have found that the IMF may be well 
parameterised by a single power-law (iSalpeter, ,1955.) , a multi- 
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segment powerlaw (e.g. iKroupa et alj|l993l : lKroup3l200l here- 
after, "Kroupa-like" IMFs), or a lognormal (e.g. iMiller & Scald 
1 19791 : Ichabrieill2003h . Generally, observations have found at most 
mild variations in the IMF within the Galaxy tBastian et alji201(]i) . 



From a theoretical standpoint, a broad range of models for the 
origin of the IMF exist. The central ideas behind the se range from 
giant molecular cloud (GMC) Jeans mass arguments jLarsonl 19981, 
EOO?; 'Bate & Bonnell' '2005'; 'T umlinsonI |2007| : lElmegreen eT^ 
^008; Narayanan & Dave 201 2,) to p rotostellar feedbac k mod- 
els liSilklll995l: I Adams & Fatuzzol [l996: Bate 2009; Krumhojd 



201 ih to models of a turbulent ISM (.Padoan & Nordlund .200^; 
Hen nebelle & Chabrieil |2008| ; lHoDkinsfl2o'l2h . While the driving 
physics behind the origin of the IMF varies substantially, nearly 
all theories suggest that the IMF should var y with changing physi- 
cal conditions, e ven if relatively weakly (c.f. lElmegreen et al.l2008l ; 
lKrumholj|201ll) . Though there is currently no unambiguous direct 
evidence for s ystematic variation s of the IMF with environmen- 
tal conditions jBastian et al ] l2O10h . the last decade of observations 
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have provided tantalizing indirect evidence for two emerging trends 
in IMF variations. 

Tiie first suggests tfiat tfie IMF may be bottom-heavy (de- 
fined as an excess of the ratio of low mass stars to high mass 
stars with respect to a Milky Way IMF) in present-epoch early- 
type galaxies. Complementary observational techniques have ar- 
rived at this conclusion, adding to its robustness. For example, re- 
cent observations of gravity-sensitive stellar absorption lines (such 
as FeH [the "Wing-Ford" band], Ca II and Na I among others) com- 
bined with new empirical stellar libraries of cool stars jRavner et alj 
I2OO9I) have been able to distinguish K and M dwarfs from K 
and M giants, and f ound the IMF to be bottom-heavy in z ^ 
early-type galaxies ( Wing &. Ford 1969; van Dok kum & Corrr^ 



20ldl2012|-|conroy & van Dokkum.2012a,b : Spiniello et alj20i; 



Ferreras et al l l2012i ; ISmith et alj|2012l)~ Meanwhile, observed con- 
straints on the stellar mass to light ratio from stellar kinematics have 
arrived at the same conclusion (e.g . Auger et al. 2010; Treu et al. 
2010';'Spiniello et al."'201 l';'Cappellari et al. 20 12iJbl ; lBrewer et alj 
2012 ; .Dutton et al., ,2012 ; Tortora et al, 2012) . Interestingly, both 



methods find that the IMF varies systematically in a manner that 
more massive galaxies have increasingly bottom-heavy IMFs. 

If elliptical galaxies are the end-product of massive starbursts, 
then the implication is that the IMF may show an excess of low- 
mass stars in heavily star-forming environments. This is, however, 
in direct contrast to a number of indirect measurements of the IMF 
in starbursts which all argue that the IMF may be top-heavy in these 
environments, rather than bottom-heavy. Evidence for this second 
trend in IMF variations comes from a variety of complementary 
measurements. 

For example, iTacconi et iD ( l2008h performed simultaneous 
modeling of the stellar masses, CO-H2 conversion factors and stel- 
lar initial mass functions of z --^ 2 Submillimetre-selected galaxies 
(SMGs), the most rapidly star-forming galaxies in the Universe, 
and found a best fit model with a mass to ligh t ratio roughly half 
that of a standard Kroupa IMF. Similarly, IvanDokkum (2 00a) sug- 
gested that the colour evolution and mass to light ratios for z ~ 1 
early-type galaxies may best be explained by a bottom-light IMF. 
More indirect evidence comes from models that have appealed to a 
top-heavy IMF at z ~ 2 to resolve the noted discrepancy between 
the observed SFR-i\f, relation and simulated one (Dave 2008). 

Some evidence for a top-heavy IMF in h eavily star-forming 
sys tems at z ~ exists as well. For example, lRiekeetal.Nl993l) 
and lForster Schreiber et ai]( l2003h suggest that the I MF may have a 
turnover mass a factor ~ 2— 6 larger than traditional Kroupa ('2002') 
IMF in the nearby starburst galaxy M82 (though see Satyapal et al. 
Il997h . Similarly, a simultaneous fit to the present day AT-band 
luminosity density, observed cosmic star formation rate density 
and cosmic background radiation by iFardal et alj ( l2007h suggest 
a "paunchy" IMF wherein there is an excess of intermediate mass 
stars. Even in the Galact ic Centre, observations sugge s t that a top- 
heavy IMF may apply jNavakshin & SunvaevI |2005| ; IStolte et al] 
I2OO5I) . While all of these observations can be explained without 
the need for IMF variation|^ it is compelling that when a variable 
IMF is inferred for a high-star formation rate surface density envi- 
ronment, it invariably tends toward top-heavy. 

In short, while the evidence is by no means firm, deviations 
from a Milky Way IMF appear to follow two trends: in present- 
day ellipticals, the IMF tends toward an excess of low-mass stars. 
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whereas in heavily star-forming galaxies from z ~ 0-2, the IMF 
appears to have an excess of high-mass stars. Simultaneously un- 
derstanding the origin of both of these trends presents a challenge 
for any theory of the IMF. 

One possibility is that systematic variations in the IMF owe 
to varying physical condition s in the star-forming molecular ISM. 
Jeans ( 190i) and others (e.g.lLarsonli2 005; Bate & Bonnell 200j; 
iKlessen et I200I iTumlinsonI \im\ INaravanan & Pavel 120121) 
have postulated that the characteristic stellar mass formed may re- 
late to the fragmentation properties of the parent molecular cloud 
(this is commonly known as the "Jeans Conjecture"). In this pic- 
ture, the temperature and density of the molecular cloud {Mj ~ 
y3/2^^i/2^ determine the minimum fragmentation mass of the 
cloud, and thus set the stellar mass scale. 

In this paper, we explore the cosmic evolution of the IMF in 
galaxies under the Jeans Conjecture to examine whether a stellar 
IMF that varies with the Jeans masses of molecular clouds can ex- 
plain the tentative evidence for bottom-heavy IMFs in early-type 
galaxies as well as the purported top-heavy IMF in high-star for- 
mation rate density galaxies. It is important to note that we neither 
argue for or against the aforementioned observed IMF trends, but 
rather take them as a given. Similarly, we make no claims regard- 
ing the validity of the Jeans Conjecture, but rather seek to assess 
its plausibility by examining the cosmological consequences of a 
stellar IMF that varies with the Jeans mass in molecular clouds. In 
§|2l we describe our methods, detailing our galaxy formation mod- 
els as well as model for the molecular ISM. In § [3] we report our 
main results, and in §|4l we provide discussion. Finally, in §[5] we 
summarise. 



2 GALAXY EVOLUTION AND ISM MODELS 
2.1 Summary of Methods 

Our principle goal in this work is to model the cosmic evolution of 
the IMF under the assumption that the IMF varies with the Jeans 
properties of molecular clouds. Because resolving molecular cloud 
scales (~ 10s of pc) in c osmological volumes is intractable for all 
but a few galaxies (e.g. Christe nsen et al. we have devel- 

oped a multi-scale methodology in order to model the mean ISM 
properties for an ensemble of galaxies. 

Our first goal is to parameterise the mean molecular ISM 
physical properties in terms of a global property. To do this, we 
employ a large suite of relatively high-resolution smoothed parti- 
cle hydrodynamic (SPH) simulations of isolated disc galaxies and 
galaxy mergers in evolution. These galaxies sample a range of 
galaxy masses and merger mass ratios, with the main goal being 
to simulate a large dynamic range of physical conditions. These 
models can marginally resolve the surfaces of the most massive 
giant molecular clouds (~ 70 pc), and employ a fully radiative 
sub-resolution model for the physical and chemical state of GMCs 
below the resolution limit. We consider the main physical processes 
that determine the thermal structure of the giant molecular clouds, 
and thus are able to resolve their Jeans properties. As we will show, 
the high-resolution SPH galaxy evolution simulations show that 
the average Jeans mass in a molecular cloud scales well with the 
galaxy's star formation rate. 

We then aim to understand the cosmic evolution of the Jeans 
properties of GMCs, and hence the stellar IMF. We model the cos- 
mic evolution of the physical pr operties of galaxie s utilising the an- 
alytic methodology described in Dave et alj ( l20l3) . This methodol- 
ogy describes the evolution of the galaxy star formation rate, stellar 
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mass growth and metal enrichment while assuming that galaxies 
live in an equilibrium between gas inflows, outflows and consump- 
tion via star formation. The results of the analytic models provide 
a good match to the evolution of the physical properties of bona 
fide hydrodynamic cosmological simulations, but at a substantially 
reduced computational cost. The main output from these models is 
the star formation histories for galaxies of various masses. Com- 
bined with a model for how the average Jeans properties of molec- 
ular clouds vary with the galaxy-wide star formation rate (as in- 
formed from the high-resolution hydrodynamic galaxy evolution 
simulations), these models then describe the cosmic evolution of 
the IMF under the Jeans Conjecture. 

Finally, we combine these results with population synthe sis 
models (FSPs dConrov et al.ll2009ll201(]| ; IConrov & Gunn|[2oTol) in 
order to explore the observable properties of stellar populations un- 
der a varying IMF. We utilise these to compare directly against ob- 
servations. 



is modeled as multi-phase, with cold, st ar-forming clouds embed- 
ded in a hotter, pressure-con fining phase jMcKee & OstrikeJ 19771 ; 
ISpringel & Hernauistll2003 l). These phases can exchange mass via 
supernova heating of cold clouds and radiative co oling of the hot- 
ter phase. Stars form according t o a volumetric jSchmidtl ( Il959h 
power-law relation with index 1.5 ( lKennicutllll998iJ) . The normal- 
isation i£_setJomatch_the local surface density Esfr — Egaa re- 
lation l lKennicutt|[T998 bV We note that while both the slope and 
normalisation of the Kennicutt-Schmidt sta r formation law are the 
subject of heavy debate (see e. g. Krumh olz & ThompsonI 1 2007 
iBigi el et alJl2008 l:'Naravanan et al. 2008^1^" 
Ostriker et alj I2O IO; Ostriker & ShettJ 



20111: 



I20 Hal : [Narayanan & Daveil2012i : Kennicutt & Evans 2012, for just 



a few examples), we showed in lNaravanan et al.. (201 lb.) that the 
physical conditions in the ISM are relatively insensitive to these pa- 
rameters so long as power-law indices between ~ 1 — 2 are adopted. 



2.2 Hydrodynamic Galaxy Evolution Models 



We simulate the hydrodynamic evolution of a large li- 
brary of idealised isolated disc galaxies and galaxy merg- 
ers using a modified version o f the p u blicly availabl e code 
GADGET- flSnringel & Hemguistl 120031 : ISpringel et aj] l2005l : 
ISpringell UmSth The galaxies range in baryonic mass from ~ 
IO^^Mq to ~ 1O ^^M0, and their physica l properties are listed in 
full in Table Al of lNaravanan et all j2012bl) . We simulate mergers 
between these discs with mass ratio 1 : 10, 1:3 and 1:1. Sim i larly, th e 
initial orbital configurations are given in lNaravanan et al] ( 1201 2bl) . 

The disc galaxies are exponential and initialised according 
to the IMo et al] j 19981) mode l . They are embedded in a live dark 
matter halo with a iHem guisti fT99of) density profile. Mergers sim- 
ply involve two discs initialised as such. The halo concentration 
and virial radius of the galaxies are motivated by cosmological A^- 
body simulations, and are redshift dependent IBullock et al.ll200l] : 
iRobertson etai]|2006l) . We simulate galaxies scaled for z=0 and 
z=3. The gravitational softening lengths are ~ 100 h"^pc for 
baryons, and ~ 200 h^^pc for dark matter 

The z=0 galaxies are initialised with a 40% baryonic gas 
fraction, and z=3 galaxies with an 80% baryonic gas fraction. 
These gas fraction s are potentially large r than those observed (e.g . 
iDaddi et alj [2010l : iTacconi et al.Tl2010l : iNaravanan et alj l2012al) . 
However, because we do not include any form of gas replenish- 
ment from the intergalactic mediu m (IGM) in these idealised sim- 
ulations (e.g. iMoster et alj 1201 lal lbl). we require large initial gas 
fractions in order to have substantial lifetimes of moderately gas- 
rich galaxies during the evolution of the merg er. When galax- 
ies merge, this is necessary to driv e a starburst dNaravanan et alj 
l20I0bl : iHavwardetallbOIll l2012al) . In any case, the exact choice 
of initial gas fraction is not terribly relevant as we principally sim- 
ulate the evolution of idealised galaxies in order to build a sample 
of galaxies with a large dynamic range of physical conditions. The 
gas is initialised as primordial. A mass fraction of stars (utilising a 
Salpeter IMF) are assumed to die instantly upon each star forma- 
tion event, and enrich the surrounding ISM with metals assuming 
instan taneous recycling and a yield of 0.02 (Springel & Hernquislj 
l2003h . 

Within the hydrodynamic evolution of the galaxies, the ISM 
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2.3 Molecular Gas ISM Specification 

We determine the physical and chemical properties of the molecu- 
lar gas in the SPH galaxy evolution simu lations in post-processing . 
We describe the methodology in full in INaravanan et alj ( 12011 
though summarise the relevant details here. 



2.3. 1 Cloud Chemical and Physical State 

We first project the physical properties of the SPH particles onto 
an adaptive mesh using the SPH smoothing kernel. The base mesh 
is 5^ spanning 200 kpc on a side, and recursively refine in an oct- 
subdivision based on the refinement criteria that the relative density 
of metals should be less than 0. 1, and that the l/-band optical depth 
across a cell be less than unity. The smallest cells refined to ~ 70 
pc across, just resolving massive GMCs. 

GMCs within the cells are modeled as spherical and isother- 
mal. The H2 fraction is determined by balancing photodissociation 
rates by Lyman- Werner band photons against H2 formation rates 
on dust grains. For simplicity, we utilise the analytic prescription of 
iKrumholz et"ai] l l2008ll2009a|jbh which assumes equilibrium chem- 
istry. A comparison of this model against full radiative transfer 
models suggests that this approximation is reasonable for metal- 
licities above Z > O.OIZq, well below the metallicities of any 
galaxy modeled here. 

The GMCs are assumed to be of constant density. In order 
to treat GMCs that reside in large cells in the adaptive mesh (and 
avoid unphysical conditions), we establish a floor surface density 
of 100 Mq pc~^. This value is motivated by observations of Lo - 
cal Group GMCs dBolatto et al1l2008l : iFukui & Kawamurdl201(j|) . 
Generally, in heavily-star forming systems, the mass-weighted sur- 
face density of G MCs is much larger than this, and the floor value 
is never reached dNaravanan et alj I2OI Ibl) . For these "resolved" 
GMCs, the cloud is assumed to occupy the entire cubic volume 
of the cell that it resides in. In lower SFR systems, however, the 
bulk of the GMCs have Eh2 ~ 100 Mq pc~2. 

With the H2 mass and cloud surface density established, 
and an assumption regarding the GMC geometry (either spher- 
ical[cubic] if the cloud is unresolved[resolved]), the density 
is known. GMCs are known to be supersonically turbulent 
dMcKee & Ostrikeill2007l and references therein.). We allow for 
the turbulent compression of gas by scaling the volumetric densi- 
ties of the GMCs by a factor e^f^"^ , where numerical simulations 
show that for supersonic turbulence: 
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ln(l + 3MfD/4) 



(1) 



where Mip is the one-dimension al Mach num ber ([Ostriker et al 



I2OO1I: IPadoan & Nordlundl llOOl though see iLemaster & Stone 
(I2OO8I)). As we will show, the temperature calculation of the GMCs 
is dependent on the densities. Because solving for the density com- 
pressions and temperatures simultaneously for many millions of 
cells across thousands of galaxy snapshot realisations is a compu- 
tationally prohibitive task, we assume the temperature of the GMC 
is 10 K for the sound-speed calculation. 



2.3.2 Cloud Thermal State 

The other pertinent variable for calculating the Jeans properties of 
the model GMCs is the gas kin e tic tem pera ture. Our model i s base d 
on that laid out byl Goldsmifl] i200lh and lKrumholz et al.l j201lh . 
The dominant gas heating terms are by grain photoelectric effect 
at a rate per H nucleus Fpc, cosmic ray heating at a rate Fcr, and 
cooling via either Cn fine structure line cooling, or ^^CO rotational 
line cooling at a rate Aune- Finally, there is an energy exchange 
between the dust and gas at a rate We assume thermal balance 
between the gas and dust which gives the following equations that 
are solved simultaneously by iterating on the temperatures of the 
gas and dust: 



F„ 



Fcr — Aiin 



Fdust — Adust — ^' 







(2) 
(3) 



We assume that the photoelectric heating rate is attenuated by 
half the mean extinction of the cloud, such that: 



Fpe = 4 X 10" 



Goe 



-]VH<Td/2 



ergs 



(4) 



where Gq is the FUV intensity relative to the Solar neighborhood, 
and (Jd is the dust cross section per H atom to UV photons. For 
simplicity, we assume that the Gq ~ 1 and (id = 1 x 10^^^ cm~^, 
though we note that tests that involve scaling the Habing field by the 
SFR make little differenc e on the thermal properti es of our clouds 
owing to strong shielding l lNaravanan et al.|[201 Ibh . 
The cosmic ray heating rate is given by: 



Fcr = C'lJCR s ^ 



(5) 



where C,' is the cosmic ray ionisation rate (here, assumed to 
be 2 X 10~^^Z's~^), and gcR is the thermal energy increase 
per cosmic ray ionisation. For H2, gcR ~ 12.25 eV (though 
note that thi s value is quite unce rtain; see discussion in Ap- 

f iendix A4 of [K rumholz et aljllollh . and for HI, gcR ~ 6.5 eV 
Dalgamo & M cCrav 19?!). 

An important aspect of these calculations is that we assume 
the cosmic ray heating rate scales linearly with the galaxy-wide 
SFR, anchored by an assumed Milky Way SFR of 2 MQyr~^ 
jRobitaille & Whitnevll2010l ; IChomiuk & Povichll201 iV While the 
exact scaling between the cosmic ray ionisation rate and SFR is un- 
known, some observational evidence exists suggesting a linear re- 
lationship. Observations of M82 by the VERITAS group, as well as 
7-ray observation^ of the Galaxy, Large Magallenic Cloud, NGC 
253 and M82 by the FERMI group all suggest a linear relation be- 
tween the galaxy -wide SFR and cosmic rays jAcciari et al]|2009l : 
lAbdoetalfcoiob . 

The dust temperature is calculated via full 3D Monte Carlo 



7-rays can arise as the result of cosmic ray interactions with H2, via 



radiative radiative transfer. We utilise StJNRISE, a publicly avail- 
able dust radiative transfer simulation package {jonssoni 1 20061 : 
Ijonsson et ai]|2006l . I2OIOI : LFonsson & Primack"201oV The detailed 
algorithms are described in[Jonsson et al. (2010), and we only sum- 
marise here. 

The sources of light are stellar clusters and an AGN which 
acts in the hydrodynamic simulations as a sink particle that 
accretes according to a Bondi-Hoyle-Lyttleton parameterisation 
Bondi & Hovlelll944l). The stars emit a STARBU RST99 spectrum 



Leitherer et alJI 1 9991 : 1 Vazquez & Leithereill2005t) . where the ages 
and metallicities derive from the SPH calculatio ns. The AGN 
emits a luminosity-dependent SED based on the iHopkins et alj 
( l2007h templates that derive from observations of unreddened type 
I quasars. In practi ce, the AGN does not im pact our results signifi- 
cantly (see tests bv lNaravanan et alj20I Ibh . The radiation traverses 
the ISM, and is scattered, absorbed and re-emitted until the dust 
temperatures converge. The dust and radiation field are assumed 
to be in radiative equilibrium. The dust mass is set by assuming a 
constant dust to metals ratio of 0. 4, comparable to th at of the Milky 
Way (Dwek 1998; Vladilo 1998; ICalura et alj2008l) . and we utihse 
the Weingartner & Draine ( 2001) i? = 3.15 dust grain models with 
Ibraine & Li (20_Q7) updates. 

For disc galaxies at low-redshift, we assume that young (< 10 
Myr) stellar clusters are embedded in photodissociation regions and 
HII regions for 2 — 3 Myr. In this case, the STARBURST99 spectrum 
is replaced by SEDs derived from MAPPI NGSIII photoionisation 
calculations dCroves et al.l200^ ; |jonsson et al..2010,). The covering 
timesc ale is motivated by parameter- space surveys bv ljonsson et alj 
( l2010h who showed that this parameter choice produces synthetic 
SEDs of 2=0 disc galaxies comparable to those observed in the 
Spitz er Infrared Nearby Galaxy Survey (SINGS; iKennicutt et alj 
l2003h . 

For galaxy mergers at 2=0 and all galaxies at 2=3, the stel- 
lar densities become so high that stellar clusters overlap. In this 
case, using MAPPINGSIII no longer makes sense as the super- 
stellar clusters (sometime as massive as ~ 10* ''^0) saturate the 
MAPPINGSIII photoionisation models. In these cases, we assume 
the cold ISM is a uniform medium with a volume filling factor of 
unity. Observations of nearby starburst galaxies may support this 
picture for the cold ISM structure in heavily star- forming galaxies 
dDownes & SolomoiJll998l ; ISakamoto et al.lll999h. Bey ond this, in 
tests run in the Appendix of Narayanan et alj 1 201 Ibi) . we found 
that the thermal and physical properties of the molecular ISM were 
generally robust against the choices made for the PDR clearing time 
or ISM volume filling factor. 

The line cooling is assumed to happen either via CII or CO. 
The fraction of hydrogen where carbon is mostly in the form of CO 
is we ll approximated by CWolfire et al..,20ICl : .Glover & Mac LowJ 
l20Tlh : 



/co = /h2 X e 



-4(0. 53-0. 0451n Q_-__0.0971nZ')/Av 

rejj/cm 



(6) 



pion-decay. 



When this fraction is greater than 0.5, we assume the cooling hap- 
pens via CO rotational lines; otherwise, the cooling is dominated 
by CII emission. We assume a constant carbon to H2 abundance of 
1 X 10"* X Z' (Lee et al. 1996), where Z' is the metallicity with 
respect to solar. 

To determine the cooling rates, we utilise one-dimension al es- 
cape probability calculations dKrumholz & ThompsoiJ [2007b . The 
level populations of the molecule are assumed to be in statistical 



Cosmic Evolution of IMF and Bottom-Heavy Ellipticals 5 



equilibrium, and determined tlirough the rate equations: 



fu 



(7) 
(8) 



where C are the coUisional ratefl A^i are the Einstein coefficients 
for spontaneous emission, / the fractional level populations, and 
Pui is the escape probability for transition u ^ I. The rate equa- 
tions can be rearranged as an eigenvalue problem, and solved ac- 
cordingly. 

The escape probabil ity, /?ui can b e approximately related to 
the line optical depth Tui l lDrainellioilh via: 

1 + 0.5r„i 
and the optical depth is: 



gi 16(27r)3/2cr \ fig„ 



(10) 



where Q is the abundance of CO with respect to H2, gi and gu are 
the statistical weights of the levels, A''h2 is the column density of 
H2 through the cloud, X^i is the wavelength of the transition, and 
a is the velocity dispersion in the cloud. 

The ID velocity dispersion, a is calculated by: 



a = max(crcoii,o-v 



(11) 



Here, ctccii is the mean square sum of the subgrid turbulent ve- 
locity and the resolved non-thermal velocity dispersion. The tur- 
bulent ve locity is calculated fro m the external pressure from the 
hot ISM teobertson et alj|2004[) . though we impose a ceiling of 
10 km s~^ which is informe d by simulations of turbulent en- 
ergy driving and dissip ation jPibetalJllOO^ : Ijoung etaljlioosl : 
lOstriker & Shettvll2oT l|). The resolved non-thermal component is 
calculated as the mean of the standard deviation of the velocities 
of the nearest neighbour cells in the three cardinal directions. A 
floor velocity dispersion of a^ii is established (typically for poorly 
resolved GMCs) by assuming the GMC is in virial balance with 

Qvir = 1, with avir = 5(7^;^ i?/ (GM) . 

With these elements in place, we then iterate Equations [2] and 
|7]-[l0]in order to determine the equilibrium gas temperature. 

Finally, we note that our model does not include any effects of 
X-ray heating of gas from either massive st ars or accreting compact 
objects (e.g. Hocuk & Spaans|[2oTol.l201lh . This may increase gas 
heating beyond the contribution of gas-dust coupling, or cosmic 
rays in rather extreme environments. 

2.4 Cosmological Models 

As will be discussed in §(3) the temperature (and more generally, the 
average Jeans mass of a GMC) in a galaxy scale with the galaxy- 
wide star formation rate. In order to understand the cosmic evolu- 
tion of the physical properties of galaxies, and therefore the stellar 
IMF , we must turn to c osmological galaxy evolution models. 

lDaveet"aLl ( l2012h presented analytic models that describe the 
cosmic evolution of the stellar, gaseous and metal content of galax- 
ies under the condition that they lie in a slowly evolving (with 



* The rat e coefficients are ta ken from the Leiden Atomic and Molecular 
Database jSchoier et alj2005l) . 



redshift) equilibrium state between gas infall from the IGM, gas 
outflows, and gas consumption from star formation. These models 
compare quite well to the results from cosmological hydrodynamic 
simulations (e.g. Oppenheimer et al. 2010), though have the practi- 
cal advantage of being computationally inexpensive. We therefore 
utilise these and summarise the relevant aspects of the model here. 

In this model, the star formation rate is determined by the 
gas infall rate minus the outflow rate, which is given by the mass 
loading factor t) times the SFR. The gas infall rate is mitigated 
by a preventive feedback parameter C, < 1, while it is aug- 
men ted by recyc ling of winds governed by the recycling parameter 
az jpinlator & D ave 2008), which is the ratio of the infalling gas 
metallicity to the ISM gas-phase metallicity. The star formation rate 
is given by 



SFR = 



(l + 77)(l-az)' 



(12) 



Hence there are three parameters in this model, rj, (, and az, and 
their depen dences on mass (u sually halo) and redshift. We take 
Mgrav from lDekel etall ( |2009|) . 

Our choices for these three baryon cycling parameters are 
guided by numerical simulations and observational results. We as- 
sume rj oc Mj^jjj^^, as advocated bv ISehroozi et alj l l20I2ah and, 
as we show in Figure [T] results in a model stellar mass to halo 
mass ratio comparable to what is inferred from abundance match- 
ing. W e take az from a cr ude parameterization based on simula- 
tions of iDave et al] ( l20Ilah 



i(Afhaio/lO''Af0,l)) 



2/3 



(13) 



These analytic model choices approximately reproduce simulation 
results for key galaxy evolutionary properties such as the star for- 
mation history, metallicities, and gas fractions. 

For this work, a key issue is quenching, because as we will 
show, low-level star formation in massive galaxies after quench- 
ing is an important element of our model. We include a number of 
sources of preventive feedback in the model (Q. 



C — Cphoto X Cgrav X Cwinds X Cqucnch 



(14) 



where ("photo represents heating of cold infalling gas by 
photoionising radiation, and is principally important at low 
masses. Here, we use the par a meter isation advocated by 
iGnedinI ( l2000h and lOkamoto et al] ( |2008|) . such that Cphoto oc 
[l + l/3{M/M^)-^]~^'^ with given in iDave et al] tOlj. 
Figure I). Cgrav parameterises the suppression of gas inflow ow- 
ing to virial shocks, and is p roportional to (1 -I- z)"'^* x M^^^^^ 
jFaucher-Giguere et alj20I li) : while the simulations that determine 
this parameterisation for Cgrav d o not include metal line cooling, 
results with metal cooling from iDave et all ( l20Ilbl) confirm this 
scaling, (winds represents the impact of energy input by winds, 
and scales as ~ e' 

V^fhaio dOppenheimer et al.ll20icl) . Generally 
speaking, winds and photosuppression of infall do not impact the 
massive galaxies that concern us in this study. 

Cqucnch represents putative quenching mechanisms (e.g. feed- 
back from an active galactic nucleus [AGN]) not encompassed in 
the other forms of feedback. For massive galaxy evolution, (quench 
tends to have the biggest impact on the SFH which can impact the 
exact form of the stellar IMF at 2 ~ in our model. Here, we 
investigate our results in t erms of two poten tial quenching mod- 
els. For the first, we follow lOave et all l l2012h . and model (quench 
with a similar functional form as (photo, but with a characteristic 
mass of A/q — 10^^'^Mq. The second model appeals to recent 
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Figure 1. M,-Mhaio relation for model galaxies with the two quenching 
models employed in this work. The red and blue lines show the two dif- 
ferent quenching models employed in this work, while the yellow band 
shows observational constraints presented bv lBehroozi et all )2012bl l3). The 
left panel shows z ~ galaxies, while the right shows galaxies at 2 ~ 2. 
The gr een crosses in the right panel show observational constraints by 
IWake et al.1 ^2011 ). At low-redshift, galaxies that form in our equilibiium 
model show excellent coiTespondence with observed galaxies. Beyond this, 
the different quenching models produce similar galaxies, owing to the fact 
that the bulk of stars are formed prior to quenching. At higher-redshifts, the 
models may be discrepant from observations at low mass, though we note 
that observations of stellar mass functions at z ~ 2 are still relatively un- 
certain. In any case, the correspondence for massive galaxies (the galaxies 
that concern us in this work) is quite reasonable, even at high-2. 



15 



work by iBehroozi et al ] l l2012ah which utilises observational con- 
straints on the stellar masses of galaxies to model the scaling be- 
tween the cosmic star formation efficiency (cSFEfland galaxy halo 
mass. Isehroozi et al] i2012ah find that the cSFE scales as M^^^^ . 
When parameterising in terms of the stellar mass-halo mass rela- 
tion, the threshold halo mass beyond which the cSFE becomes in- 
efficient is ~ Mhaio > 10^^ Mq . B ecause Cgrav dominates over 
Cphoto and (^winds for massive galaxies, including a quenching term 
Cquench oc M^^^^'^^ for halo masses A/haio > IO^^Mq results in 



a cSFE in our models tha t scales with M^^^^J' , per the constraints 
of lBehroozi et alj j2012ah . 

The galaxies that result from the two quenching models are 
reasonably similar. In Figure [T] we plot the z=Q stellar mass-halo 
mass relation for the galaxies in both quenching models. Owing to 
the fact that the bulk of the stellar mass is built up prior to quench- 
ing, the two models results in quite similar M,-Mhaio relations. 
As we will show, while the low-level star formation history at late 
times can affect the present-day stellar IMF in our model, for plau- 
sible choices of the SFH, the results are qualitatively similar. 



^ i.e. The galaxy SFR divided by the halo accretion rate. This is to be dis- 
tinguished from either the mass of stars formed divided by the total mass 
in a molecular cloud (Af»/(M» -|- Mdoud)). or the SFR/AfH2 (i-e- an in- 
verse time scale), both of which are also commonly referred to as the star 
formation efficiency. 



This equilibrium model for the cosmological evolution of 
galaxies does not account for galaxy mergers. While merger-driven 
starbursts d o little to contribute t o the total stellar mass of a 
galaxy (e.g. iRodighiero et alj|201ll) . the added stellar content can 
affect the observed properties (e.g. the mass to light [M/L] ratio) 
if the merger is close to equal mass. Hence we explicitly include 
additional growth owing to galaxy mergers. 

We estimate the galaxy -galaxy merger rate f rom the publicly 
available model developed in lHopkins et al. Briefly, we cal- 

culate halo mass functions and merger rates at a given redshift from 
the Millenium simulation l lFakhouri et al.ll201ol) . and assign galax- 
ies to these haloes following a standard halo occupation formal- 
ism ( Conr ov & Wechsleil2009l) . After correcting for the dynamical 
friction time-delay between halo-halo and galaxy-galaxy mergers 
teovlan-Kolchin et al.ll2008l) . we then have a galaxy merger rate 
as a function of galaxy mass and redshift. The mapping between 
galaxy mass and halo mass is dependent on knowing the galaxy 
mass function as a function of redshift. For this, we utilise observed 
galaxy ma ss functions. From z = — 2, we uti l ise the mass func- 
tions from lAmouts et al.l j2007h and lllbert eFZI j2010l) as bracket- 
ing the range of observed mass functions. At 2: > 2 we utili se mass 
funct ions from Perez-Gonzalez et al. (2008) and Marchesi ni et alj 
( I2OO9I) . The principle effect of including galaxy mergers is to add 
to the uncertainty present in trends as the uncertainty in high-z 
galaxy mass functions dominates the uncertainty in these models. 
This said, as we discuss later, the trends presented in this work are 
robust, and only the normalisation tends to vary with different mass 
function choices. We have run versions of our model with each of 
the aforementioned mass functions, and discuss the range of uncer- 
tainty as we present our results. 

The combination of simple but realistic star formation histo- 
ries for individual galaxies from the equilibrium model plus the 
merger history derived from the Millenium simulation allows us 
to quickly and effectively model the evolution of the galaxy pop- 
ulation. This in turn enables us to efficiently examine the impact 
of IMF variations parameterized by SFR and metallicity on galaxy 
properties. 



3 RESULTS 



Summary of Cloud Thermal State and Scaling of IMF 
with Galaxy Physical Properties 



3.1 



Because the methods described in §|23] dictating the thermal state 
of the clouds are somewhat detailed, we find it useful to summarise 
the general trends in cloud temperatures with the intent of giving 
the reader an intuition for the different physical processes that de- 
termine the temperature in different regimes. 

Without any sources of heating, the cloud temperatures would 
cool to the microwave background temperature. The dominant 
sources of heating are gas-dust coupling when the GMC is at high 
density (n > 10^ cm^'^), and cosmic ray heating at lower densities 
(n < lO^cm"^). 

At high densities, the energy exchange between gas and dust 
becomes extremely efficient. In heavily star-forming systems, the 
gas temperature is, to first order, s et by the d ust temperature due to 
large fractions of dense gas (e.g. ljuneau et al. 2009). While cos- 
mic rays can play a significant role in high SFR galaxies (e.g. 
|PaDadoDoulosl2010l ; |PaDadoDoulos et alj201lh . in these models en- 
ergy exchange with dust plays a comparable or greater role in set- 
ting the gas temperature jNaravanan et alj2012bl) . Because the dust 
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temperature rises with the SFR (e.g. iNaravanan et alj|2010al) . the 
gas temperatures do as well. 

In low SFR systems, the mean gas density is relatively low. 
Here, cosmic rays determine the gas temperature. For a galaxy 
forming stars at a rate similar to the Galaxy (~ 2 MQyr~^), the 
resulting temper ature is ~ 10 K, comparable to observation s of 
Galactic clouds ( iBlitz et alJboOTi ; iFukui & Kawamur3l201(]|) . At 
lower SFRs, the cloud temperature can drop below this value. 

The cloud Jeans mass, which scales as T'^^^ /n^^^ therefore 
scales, on average, with the galaxy-wide star formation rate. Phys- 
ically, this dependence arises from gas-dust coupling at high den- 
sities (combined with increasing dust temperatures with increas- 
ing star formation rates), and the scaling of cosmic ray fluxes with 
SFRs at lower densities. Utilising the information from the ~ 10^ 
galaxy snapshot realisations modeled here (agai n, the entire pa- 
ramete r space surveyed is listed in Table Al of 'Narayanan et al.' 
l2012bl) . we fit the relationship between the H2 mass-weighted Jeans 
mass in our model galaxies and the galaxy SFR. The H2 mass- 
weighted Jeans mass is calculated by weighting the Jeans mass of 
each cell by the H2 mass present, and the Jeans mass is calculated 
utilising the median gas density in the cell (i.e. the density above 
which half the mass resides). The median density is calculated as- 
suming the gas has a lognormal density distribution function, with 
width given by Equation[T] This results in a relation: 



Mj 



SFR°-^M, 







(15) 



with a normalisation set at a Galactic SFR of 2 MQyr"^. 

Under the Jeans conjecture, the scaling of the average Jeans 
mass of molecular clouds has implications for the stellar IMF. We 
assume the IMF has a similar basic shape in all galaxies, and takes 
the form of a Kroupa broken power-law with a a low-mass slope of 
0.3, and high-mass slope of -1.3 (where the IMF slope specifies the 
slope of the logio(dN/dlogM)-logio(M) relation). For Milky Way 
conditions, the turnover mass, Mc, is set to 0.5 Mq, and varies 
with the SFR for other galaxies. Thus, for IMF variations, we have: 



Mc = 0.5 X 



SFR 



SFRmw 



(16) 



Under the Jeans conjecture, the IMF therefore varies with the 
conditions in the ISM. Physically, the IMF tends toward top-heavy 
in heavily star-forming galaxies due to efficient energy exchange 
between gas and dust in dense environments. Conversely, the IMF 
tends toward bottom-heavy in poorly star-forming galaxies owing 
to decreased cosmic ray fluxes. It is important to note that the shape 
of the IMF does not change in this model - only the turnover mass. 
To explicitly show this point, in Figure |2l along with a standard 
Salpeter and Kroupa IMF, we show example IMFs for a low-SFR 
bottom-heavy system (Mc = 0.25 Mq), and a high-SFR top- 
heavy system (Mc — 1 Mq). 

As we will see. Equation 1 1 61 combined with the results from 
§ nil] prepare us to analyse the cosmic evolution of the IMF under 
the assumption of the Jeans conjecture. Finally, we note that for- 
mally, the rise in Jeans masses in high-SFR galaxies can be modu- 
lated in very low metallicity environments owing to a metallicity - 
dependent dust-gas energy exchange rate jKrumholz et alj|201 ih . 
Typically, heavily star-forming systems do not remain metal- 
poor in our model for very long, however, and any metallicity- 
dependence in Equation [T6] has minimal impact on the IMF in 
massive galaxies at 2: ~ under the Jeans conjecture. Future work 
will discuss the star formation properties of low-metallicity, heavily 
star-forming systems at high-redshift. 
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Figure 2. Exa mple stellar initi al mass functions. The black solid line de- 
notes a standard fSalpeteJ jl955h IMF, whereas the other three show lKroupal 
( 2002) broken power-law forms. The blue dashed line is the typical Galac- 
tic Kroupa IMF with turnover mass Mc = 0.5 Mq- The green dash- 
dot line represents a bottom-heavy IMF (with Mc = 0.25 Mq ), while 
the red dashed line is a top-heavy IMF (with Mc = 1 Mq ). At solai' 
metallicities, the bottom-heavy IMF may represent a galaxy forming stars 
at ~ 0.2 M0yr~^, while the top-heavy IMF is characteristic of a galaxy 
forming stars at ~ 20 — 25 MQyr"'^. An important point is that the shape 
of the IMF does not change under the model assumption of the Jeans con- 
jecture: only the turnover mass. The IMFs ai'e normalised arbitrarily by the 
integral J M(dN/dlogioM) dM. 



3.2 The Meandering of the IMF over Cosmic Time 

In Figure [3] we plot the cosmic evolution of the halo mass, M*, 
SFR, IMF characteristic mass (c.f. Figure |2]l, and K-band mass 
to light ratio (normalised by what one would see for a standard 
Kroupa IMF) for three model galaxies for both quenching models. 
The model galaxies are chosen to have the same final stellar mass. 
We defer discussion of the mass to light ratio evolution until §[33] 
The characteristic star formation history of galaxies is, to first 
order, determined by the gas accretion rate from the IGM which 
is dictated by gravity. Fits to models from various groups all find 
accretion rates dependent on both the halo ma ss, as wel l as red- 
shift via Mgray PC M°-°^~°-^^ X z^-^-^ (e.g. iDekel et ; 
iFakhouri et ani201(]| : lFaucher-Giguere et alfcOlC 



rai]|2009l: 



In this picture, more massive galaxies on average have hig her 
peak SFRs, and peak at earlier times (e.g. iNoeske et alj|2007al lbl). 
As an example, while a (^=0) Milky Way sized-halo peaks in SFR 
between z — — 1, and has a peak SFR of ~ 10 — 20 MQyr~^, 
a galaxy in a (z = 0) ~ lO" Mq halo peaks at ~ 200 - 300 
MQyr~^ at z ~ 2 — 3. Modulo metallicity effects (which only play 
an important role for the most part at redshifts 2 > 4), the evolution 
in the turnover mass in the IMF broadly follows the star formation 
history. Turning to the third panel in Figure [3] we now highlight 
the grey horizontal line which denotes the typical Kroupa turnover 
mass of Mc = 0.5 Mq in the Milky Way. 

Massive galaxies that undergo heavy SFR periods at early 
times see an increase in their average cloud temperature^ and 



" The cloud densities of course increase as well, but the Jeans mass goes 
asMj ~T3/2/nl/2. 
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hence an increase in the characteristic turnover mass in their stel- 
lar IMFs at early t imes (e.g. top-heavy/bottom-light IMFs). In 
iNaravanan & Davi {20l5), we showed that the top-heavy IMF in 
high-z star-forming galaxies that results from a Jeans model can 
alleviate tensions associated with inferred SFRs from high-z galax- 
ies. 

As the SFRs of massive galaxies decline due to a combina- 
tion of gas exhaustion and cosmic expansion of the Universe, so do 
their cosmic ray fluxes. Below 2 ~ 2, when the cosmic microwave 
background (CMB) temperature decreases below Tomb 10 K, 
the temperatures of the GMCs may decrease below the typical ~ 10 
K value seen in Milky Way clouds. Consequently, the cloud Jeans 
masses drop and the IMF becomes bottom-heavy. For the most 
massive galaxies we consider, the transition between top-heavy and 
bottom-heavy (defined as when the turnover mass in the IMF tran- 
sitions from above 0.5 Mq to below) happens at 2 ~ 1 — 2. 

Under the Jeans conjecture, then, the IMF meanders over cos- 
mic time, generally guided by the galaxy-wide SFR. The exact 
form of the meandering is galaxy mass-dependent. For example, 
the most massive galaxies (Mhaio ~ 10^* Mq at 2 ~ 0) peak 
in their SFR, and hence their IMF characteristic mass, at 2 ~ 2-4. 
They then go through a bottom-heavy phase for roughly half a Hub- 
ble time. On the other hand, galaxies that reside in Milky Way-sized 
haloes at 2 ~ peak in their SFRs so late that they never undergo 
a bottom-heavy phase, and are only weakly top-heavy at late times 
(2 < 1). Interestingly, galaxies that are much less massive than the 
Millcy Way never have large (> 2 MQyr~^) star formation rates, 
and thus have bottom-heavy IMFs for the bulk of their lives in this 
model. 

It is important to note local peaks and valleys in the cosmic 
evolution of the star formation history from a given galaxy can 
cause deviations in the IMF from the trends presented in Figure |3] 
and discussed in this section. For example, a gas-rich major merger 
that drives an intense burst of star formation can temporarily drive 
up gas temperatures, and hence the cloud Jeans masses. This said, 
under the Jeans conjecture, on average the IMF trends will follow 
those shown in Figure|3] We now expand upon the trends discussed 
here and compare more specifically to inferred measurements of 
the IMF from observations. 



3.3 The Mass-to-Light Ratio in 2=0 Massive Ellipticals 

As discussed in the Introduction, various circumstantial arguments 
favour a bottom-light IMF in rapidly star-forming galaxies at 
high redshifts. Such argu ments are detailed more thoroughly in 
[Narayanan & Davg ( l20I2h . where we consider the implications of 
the Jeans mass conjecture on high-2 galaxies. In i]4.2l we discuss 
arguments for and against the interpretation of a bottom-light (or 
top-heavy) IMF in massive high-2 galaxies. Here, we consider the 
implications of our varying IMF for the descendents of massive 
high-2 galaxies, which today are expected to be early-type "red and 
dead" ellipticals. 

A number of works recently have found evidence for a po- 
tentially bottom-heavy IMF in passive, early-type galaxies. As dis- 
cussed in §[T] the methods for indirectly inferring the IMF include 
the usage of gravity-sensitive absorption l ine signatures combined 
with new cool star stellar libraries (e.g. IConrov & van DokkumI 
[2012b, and ref erences therein), strong lensing, and stellar kinematic 
measu res (e.g. lTreu et alj201ol ; lDutton et alj20I2l ; ICaDpellari et al] 
l20I2al lbl and references therein). While the methods are diverse, 
most results point toward a scenario in which the IMF be- 
comes increasingly bottom-heavy with increasing stellar velocity 



dispersion (e.g. 



ICappellari et al.l 



Treu et al.| I20T0I: 'Conrov & van Dokkuii] l2012bl. 
2012al : Smith et al. 2012). The inference is, of 



course, that more massive galaxies have a more bottom-heavy 
IMFs. 

At face val u e, thi s goes against our assertion in 

[Nar ayanan & Davi ( l20I2h . as well as earlier works jPavd 
2008), that progenitor massive galaxies at high 2 should have 
bottom-light IMFs. But the story isn't so simple. The stars in 
today's ellipticals are mostly old, and hence the 2 = IMF in 
ellipticals is only being probed at star masses of < 1 Mq. In 
contrast, the star formation rate indicators that point towards a 
more bottom-light IMF at high-2 trace massive stars, typically 
> 10 Mq. Therefore the bottom-heavy and bottom-light IMFs are 
only contradictory under the assumption of a single, invariant IMF 
in these systems. 

What happens in our Jeans conjecture model? Referring again 
to Figure|3] a trend in the characteristic IMF mass at 2 ~ is appar- 
ent with increasing galaxy mass. While a Milky Way mass galaxy 
peaks in its SFR around 2 ~ 1-2, and only settles to an IMF with 
il/c ~ 0.5 Mq around 2 ~ 0, more massive galaxies peak in their 
SFR at earlier times. These galaxies exhaust their fuel rapidly, are 
unable to replenish their gas owing to strong preventive feedback, 
and settle into a low SFR mode (perhaps in an episodic fashion) 
for up to ~ half a Hubble time. These galaxies produce substan- 
tial numbers of low-mass stars as they evolve. Moreover, they can 
accrete a large number of small galaxies through mergers that also 
have a more bottom-heavy IMF; for now we set this effect aside, 
but we will return to it later. 

Observationally, a key piece of evidence for a bottom-heavy 
IMF in massive early-type galaxies is the fact that the K-band 
mass-to-light ratio {M/Lk), tracing older stars, increases with 
stellar mass more quickly than expected for a standard Kroupa or 
Chabrier IMF. Therefore, we use our model to make predictions for 
A//L_fs-(o"), where a is the galaxy velocity dispersion. 

We employ stellar population synthesis models to derive the 
stellar mass to light ratio from our model galaxies. We utilise FSPS, 
a flexible, publicly avai l able stellar p opu lation synthesis code de - 
scribed in lConrov etal] j2009ll20IOl) and lConrov & GunnI ( l2010h . 
For a given galaxy, we utilise the star formation history and cosmic 
evolution of the IMF to evaluate the properties of the 2=0 stellar 
populations. We model a given galaxy by discretizing the star for- 
mation history into a series of bursts. We run FSPS assuming the 
star formation history is a single burst at each redshift, with the 
corresponding A/c derived from our Jeans model. We then sum the 
mass to light ratios, weighting by the stellar masses, of the galaxy 
over its cosmic history to determine the mass to light ratio at a given 
redshift. 

We now refer to the bottom panel of Figure|3] Here, we show 
M/ Lk normalised by what is expected from a standard Kroupa 
IMF. The masses include the contribution from stellar remnants 
(e.g. white dwarfs and neutron stars), and we discuss this further 
in § |4] As is evident, a galaxy residing in a Milky Way sized halo 
at 2 ~ has a M / Lk ratio comparable to what is expected from a 
standard Kroupa IMF. Galaxies of increasing mass go through their 
peak star formation rate epochs at earlier and earlier times, thus go- 
ing through increasingly longer bottom-heavy phases. As a result, 
more massive galaxies have larger M/ Lk ratios at 2 = 0. 

This trend is more clearly seen in Figure [4] where we exam- 
ine the 2 = mass to light ratios in our model galaxies against 
stellar velocity dispersion (serving as a proxy for galaxy mass). To 
remain consistent with the way observations are presented, we nor- 
malise M/Lk by a Milky Way ratio (here, a Kroupa IMF with 
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Figure 3. Cosmic evolution of model galaxy halo masses, M, , SFRs, IMF turnover masses, and X-band mass to light ratios for three galaxy masses. The 
left column shows the quenching model from Dave et al. (2012), while the right column shows the quenching model designed to match the constraints of 
iBehroozi etal] i20\2^) . The galaxies are picked to have the same final stellar mass for each quenching model. The thin grey line in the fourth panel denotes 
Mc = 0.5 Mq, the turnover mass for a Galactic Kroupa IM F. The mass to light ratios are normalised by the expected values for a Kroupa IMF to assist in 
comparing with observations jConrov & van Dokkuml2012bh . The small features in the AI/L ratio owe to the discretization of the SFH into individual bursts, 
as well as galaxy mergers (see text for details). The IMF for a massive galaxy may vary over cosmic time from top-heavy to bottom-heavy as follows. During 
epochs of heavy star formation, the Jeans masses of clouds increase owing to gas in clouds heated by warm dust in dense environments. These galaxies may 
exhibit top-heavy IMFs under the Jeans conjecture. Galaxies of increasing mass peak in their cosmic star formation history at earlier times, and undergo long 
periods of star formation at rates below the Milky Way's. Decreased cosmic ray fluxes during low SFR periods result in cloud Jeans masses lower than typical 
Galactic values, and hence an excess of low-mass stars in their IMFs at late times. 
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Mc = 0.5 Mq; C. Conroy, private communication). The velocity 
dispersion is calculated utilising fits to the dynamical mass ver- 
sus velocity dispersi on from Sloan Digitized Sky Survey (SDSS) 
galaxies presented in lvan de Sande et al. (201 1). 

The blue triangles and purple squares in Figure |4] denote the 
individual galaxy realisations from the cosmological simulations. 
The overall normalisation of the simulated points is not very ro- 
bust against our model assumptions. Variables such as the assumed 
Galactic SFR (which dictates when decreased cosmic ray fluxes 
will allow cloud temperatures to drop below ^ 10 K) and the as- 
sumed stellar mass function at high-redshift (which can affect the 
galaxy merger rate; c.f. i) l2.4t can slightly change the normalisation 
of the A// Lk ratio as a function of mass. The typical amount vary- 
ing these parameters within reasonable bounds changes the M/Lk 
ratios is ~0.2 dex and this region of uncertainty is highlighted by 
the grey shaded region (only noted for the Dave quenching model). 
The magnitude of the trend in M/Lk with o depends on the exact 
form the SFH (and hence the quenching model). This said, for plau- 
sible star formation histories, a general trend of increasing M/Lk 
ratio with o is evident. 

Th e points with erro r bars show data from the observational 
work of IConrov & van Do kkum (2012b). The green points are a 
sample of nearby early type galaxies (including the bulge of M3 1), 
and the red point is stacked data from 4 galaxies in the Virgo cluster. 
The data show increasing mass to light ratios with increa sing stellar 
veloci t y dispersion. Similar tr ends have b een noted_by iTreu et al] 
( l2010l) : ISmithetal .1 ( 12012) and lCappellari et al. (2 012al) as well. 

The model results are qualitatively in agreement with observa- 
tional data. They indicate that galaxies with a typical of the Milky 
Way represent a minimum in the mass-to-light ratio in our Jeans 
conjecture IMF relative to a Kroupa IMF. Galaxies at both higher 
and lower velocity dispersions tend to have increased mass-to-light 
ratios. This perhaps counter-intuitive result comes from the pro- 
tracted epoch of low star formation associated with massive galax- 
ies that populate the low-mass end of the IMF at late times, while 
in low-mass galaxies (which are typically still star-forming today) 
the low SFR's and metallicities result in a bottom-heavy IMF at all 
times. 

Quantitatively, this model cannot attain the ~ x 2 increase in 
M/Lk relative to Kroupa as seen in some galaxies, including the 
stack of Virgo ellipticals (red point). It may be possible to tweak 
model parameters to achieve this. For instance, we have assumed 
that the IMF form is always Kroupa-like; if instead the IMF also 
had a shift to a shallower slope, this would result in more massive 
stars that would disappear by today, and would increase the rela- 
tive contribution of the low-mass stars created during the low-SFR 
phases of evolution, when the IMF was more bottom-heavy. 

Galaxy mergers moderately enhance the bottom-heavy IMF in 
massive galaxies. This is because in our model, low mass galaxies 
also have a slightly bottom-heavy IMF as they never achieve signif- 
icant periods of star formation above SFRmw- As a result, mergers 
with low mass galaxies can add more low mass stars and slightly 
increase the observed M/Lk ratio. This said, we do not account 
for the increased SFR driven by galaxy mergers. Qualitatively, one 
would expect to observe a top-heavy IMF in an ongoing merger 
due to increased star formation rates and hence Jeans masses (via 
increased gas temperatures). 

A complete accounting of all these effects would require im- 
plementing a full evolutionary IMF model within a cosmologically- 
based galaxy formation model that includes the effects of mergers 
self-consistently; we leave this for future work. Nonetheless, it is 
encouraging that using reasonable models for the star formation 
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Figure 4. Stellar mass to (^C-band) light ratios as a function of galaxy 
stellar velocity dispersion for simulated galaxies (blue triangles and pur- 
ple squares), and observed gala xies (green and red circles with error bars; 
Como Y_& vanDokkurnl \2X) 1 2bl) . The blue triangles have SFHs with the 
DaveetalJi2012l) quenching model imposed, while the purple squares have 
a quenching model designed to matc h constraints on the cosm ic star forma- 
tion efficiency-Miiaio relation from iBehroozi et aD i2012al) . The M/Lk 
ratios are normalised by what is expected from a Kroupa IMF. The grey 
shaded region denotes a 0.2 dex uncertainty in the normalisation of the 
models (see text for details). Galaxies of increasing mass have their peak 
star formation rates at earlier epochs, thus undergoing long low SFR peri- 
ods. During this time, decreased cosmic ray fluxes drive down gas tempera- 
tures, and, under the Jeans conjecture, drive a bottom-heavy IMF. Galaxies 
of increasing mass have more bottom-heavy IMFs. The simulations and data 
show reasonable coiTespondence. 



history of massive ellipticals, it is possible to attain a bottom-heavy 
IMF today even though the IMF of stars forming at high-redshifts 
was more bottom-light. This embodies the key result of this pa- 
per, that bottom-heavy and bottom-light IMFs inferred in massive 
galaxies today and their progenitors at 2 ~ 2+ can be simultane- 
ously accommodated within a well-motivated hierarchical galaxy 
formation framework. 



4 DISCUSSION 

In [Narayanan & Dave! ( 1201 2h we discussed the implications of an 
IMF following the Jeans conjecture on the properties of high- 
redshift star forming galaxies. We showed that this IMF goes signif- 
icantly towards alleviating a number of difficulties in understanding 
the evolution of the SFR- A/, relation and sub-millimetre galaxies, 
and has interesting implications for the star formation law. 

In this paper we consider the implications of the Jeans conjec- 
ture on the evolution of passive, massive galaxies. While our main 
focus is their mass-to-light ratios, there are a number of other im- 
plications for massive galaxies that we now briefly consider. 
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4.1 Galaxy Star Formation and Chemical Enriciiment 
History 

For a typical massive galaxy in this model, the star formation 
history peaks at high-redshift, and once the halo mass exceeds 
the quenching mass i\/qucnch its star formation begin s to become 
strongly suppressed (see Figure 1 of iDave et alJl2012h . The phys- 
ical mechanism fo r this is purported to be AGN feedback (e.g. 
ICroton et alj|2006h , although our model is purely parameteric. The 
resulting star formation histories for different quenching models 
are shown in Figure [3] more massive galaxies have earlier and 
higher peaks of star formation, consistent with ste llar popula- 
tion data for early-type galaxies (" Thomas etZll2005h . and remi- 
niscent ofjhe_obseryationally-derived "staged" galaxy formation 
model ^Noesk e et al.ll2007al lbl). Importantly, the star formation rate 
remains finite, albeit small, even after quenching has kicked in. 
This SFH integral is important for driving top-heavy IMFs in galax- 
ies during high star formation rate epochs, and bottom-heavy IMFs 
during low star formation rate evolution up until the present-day. 

The existence of a long low-star formation rate period is 
broadly consistent with observations. Recent studies find that low 
levels of star formation and dense molecular gas is w i despread 
among z ~ early type galaxies (Young et al. 201 1 ; Xu et al.' 
l201(]|Jcrocker et aljf201ll : lDavis et al. 2011 ; Crocker et al. 2012), 
Additionally, observations by the ATLAS^'^ collaboration appear 
to tentatively suggest that ellipticals have star formation histories 
consistent with those presented here, wit h more massive galaxie s 
typically having older stellar populations dMcDermid et al.|[2012h . 



In this sense, that our model produces low-level star formation in 
early type galaxies may be reasonable, and supported by observa- 
tional evidence. 

While early-type galaxies clearly undergo some small levels 
of star formation, there is significant variance in the amount of star 
formation among galaxies, suggesting a duty cycle for this "frost- 
ing" of star formation. In our simple model, this would result in a 
less bottom-heavy AI/Lk, because the actual star formation rates 
would be higher (albeit over shorter periods of time). However, this 
would not result in a change in the trend of M/Lk unless there 
was a systematic variation in the duty cycle with a. 

Another issue relates to the chemical enrichment history of 
ellipticals, which can be significantly impacted by changes in the 
IMF. Galaxies with short star formation time scales are expected to 
have enhanced [a/Fe] ratios. This is because when star formation 
is burst-like, the enrichment is dominated by type II supemovae, 
thus the stars will be enhanced in elements that are products of the 
a process in stellar evolution. On the other hand, when star forma- 
tion occurs over an extended period, there is time for type la super- 
novae to contribute Fe to the chemical makeup of the galaxy. Ob- 
servationally, more massive galaxies appear to be more a-enhanced 
dTrager et alhood) . 

At face value, protracted periods of low star formation rate 
in massive galaxies as seen in Figure [3] might suggest decreasing 
[a/Fe] ratios with a. However, this may be compensated by the 
varying IMF. By examining chemica l evolution in g a laxy f ormation 
models that include variable IMFs, lArrigoni et alj d2010l) showed 
that more top-heavy mass functions will increase the [a]/Fe ratios 
in galaxies. Because a core feature of the Jeans model is that more 
heavily star-forming systems have increasingly top-heavy IMFs, 
this may compensate for the extended period of Fe enrichment at 
later times. Without introducing a full chemical evolution model 
(and its associated uncertainties) to these simulations, however, be- 
ing more quantitative is currently infeasible. 



In summary, observations sh ow both evidence fo r ongoing star 
formation in early-type galaxies dCrocker et alj|2011 ), older stellar 
populations on average for mor e massive galaxie s Thomas et alj 
2005; McDermid et al. 2012; P acifici etal]|2012h , and increasing 
[a/Fe] ratios in early types dTrager et al. |200(]|) . suggesting that 
the typical star formation histories calculated by our cosmological 
galaxy evolution model may not be unreasonable. However, current 
uncertainties in the detailed SFH and chemical enrichment histories 
of early-type galaxies are such that these issues cannot definitely 
confirm or rule out a Jeans conjecture origin for the bottom-heavy 
IMF in ellipticals. 

4.2 Is the IMF in Rapidly Star-forming Galaxies Top-Heavy? 

The Jeans conjecture IMF model predicts that ellipticals undergo 
an early bottom-light /top-heavy phase du ring its peak star forma- 
tion rate epoch. N arayanan & Dav^ d2012h have suggested that this 
may go some distance toward alleviating a number of difficulties 
in understanding the evolution of high-z galaxies. In the Jeans 
model, the bottom-heavy IMF in 2 ~ ellipticals arises from the 
protracted low-star formation rate period afterwards. 

There is, alternatively, a general class of models that aims 
to understand the origin of the bottom-heavy IMF in ellipticals 
by making the IMF more bottom-heav y duri ng phases of heavy 
starbursts. Bate (2009) and Krumhol^ d201ll) examined the role 
of radiative pressure in setting the characteristic mass of the IMF. 
Krumholz showed that if the characteristic mass is set by radia- 
tive feedback from young protostars, then the IMF will scale very 
weakly with pressure such that in high-pressure systems (i.e. the 
heavily star-forming galaxies that serve as the progenito rs for lo- 
cal ea rly-types) a bottom-heavy IMF results. Similarly, iHopkinsI 
d2012h showed that if the IMF is set by the sonic mass in gas, 
rather than the Jeans mass, then starbursts will show an excess of 
low-mass stars. While neither model directly examines the observ- 
able properties of a z ~ early-type galaxy population, it is likely 
that they too would show an excess of low-mass stars in accord 
with observations. In this sense, whether the IMF is top-heavy or 
bottom-heavy in starburst systems may be a discriminant in gen- 
eral classes of IMF models. It is therefore worth examining claims 
of top-heavy IMFs in high SFR systems critically. 

A number of studies have noted a mismatch between the 
observed cosmic evolution of stellar mass density and cosmic 
star formation rate density such that the inferred SFR den- 
sity prod uces a factor ~ 2 more stars at z ~ than are 
measured feopkins & Beaconill2006l ; IWilkins et aillloOS) . While 
iNaravanan & Dav^ d2012l) showed that an IMF that varies with the 
Jeans masses in galaxies could go some distance toward recon- 
ciling these differences (producing only a mismatch in the stellar 
masses at z ~ of a factor 1.3 instead of 2), some observationa l 
studies either fail to find such a mismatch (e.g. ISobral et alj|2012h . 
or propose non-IMF based solutions to the apparent discrepancy 
dReddv & Stei"dell2009l ; IStark et alj|2012l) . 

Another commonly invoked candidate for a top-heavy IMF 
in starburst galaxies is the most heavily star-forming galaxy pop- 
ulation in the Universe, the z ~ 2 Submillimetre-selected popula- 
tion. The arguments are primarily theoretical. Simply, most stan- 
dard galaxy formation models have been unable to account for the 
full observed distribution of SMGs without resorting t o a top-heavy 
IMF. The most extreme example was reported by iBaugh et alj 
(f2O05) who suggested a flat IMF may be necessary to reconcile 
the observed SMG number counts, Lyman-break galaxy popula- 
tion and present-day stellar mass function. However, recent nu- 
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merical simulations by iHavward et alj ( 1201 Oh and iHavward et alj 
( l2012bh show that if one accounts for a combination of effects, 
including small scale dust obscuration in highly-resolved simula- 
tions, the contribution of isolated galaxies, and the contribution of 
galaxy pairs, the observed SMG number counts may be accounted 
for even under the assumption of a standard Kroupa IMF. 



Ivan Dokkud ( [200^ utilised the fact that the luminosity and 
colours of a galaxy may evolve differently, depending on the IMF, 
and found that a bottom-light/top-heavy IMF best fit the evolu- 
tion of the M/Lb and U — V colours for < z < 1 early- 
type galaxies in clusters. However, when accounting for potential 
variations in population synthesis models, frosting of young stellar 
populations in the gal axies ( , T r ager et a l. 20 08), and the structural 
evolut ion of galaxies jHolden et alfcoioh . lvan Dokkum & Conrovl 
( I2OI2I) found that the same observations could be consistent with a 
Salpeter IMF. 

In general, the purported IMF variations toward top-heavy are 
all fairly mild, and while it is tempting to solve a number of appar- 
ent problems at high-z with a mild shift in the IMF towards bottom- 
light/top-heavy, there could be other explanations for each of these 
observations. This said, it is interesting that all inferred IMF vari- 
ations in heavily star-forming systems go in the direction of being 
more top-heavy/bottom-light. Whether heavily star-forming galax- 
ies have a top-heavy or bottom-heavy IMF should be a direct dis- 
criminant between general classes of IMF theory. 

Whether present-day ellipticals have gone through an early 
bottom-light phase may be distinguishable via their measured stel- 
lar mass to light ratios. Both bottom-heavy and bottom-light IMFs 
can cause increased M/ Lk ratios from the average Galactic value. 
The former increases M/ L via increased numbers of low-mass 
stars, while the latter via a contribution from stellar remnants (i.e. 
white dwarfs and neutron stars). The M/ Lk ratios from FSPS that 
we presented account for both of these effects. 

In principle, one could use remnants to potentially distinguish 
whether the IMF for massive galaxies has gone through a top-heavy 
phase at early times. If galaxies undergo a top-heavy phase, then 
their dynamical masses will have a larger contribution from stel- 
lar remnants. On the other hand, if the bulk of the stellar mass is 
built up in a bottom-heavy phase, then the contribution from rem- 
nants to the dynamical mass will be minimal. For a Chabrier IMF, 
the current mass in remnants is only a modest fraction of the living 
stellar mass, and given the relatively mild variations in the IMF pre- 
dicted here, it is unlikely that remnants would dominate the mass 
even in our most extreme cases. Nonetheless, offsets between direct 
M/ L measurements from stellar spectral features (which measure 
the masses of living stars) and dynamical modeling that includes 
the effect of remnants may potentially provide a constraint on this 
model. 

The main idea of this paper is that a bottom-heavy IMF in to- 
day's ellipticals is consistent with a top-heavy IMF at earlier times 
in heavily star-forming systems under a Jeans model for the IMF. 
This said, a bottom-heavy IMF at all epochs would obviously more 
straightforwardly result in a bottom-heavy IMF in today's ellip- 
ticals, at the cost of potentially exacerbating possible issues with 
high-redshift galaxies. Observations that distinguish whether the 
IMF varies in he avily star-form ing systems, or as a function of for- 
mation age (e.g IZaritskv et ali .20I2) may prove to be a valuable 
model discriminant. 



4.3 The Validity of the Jeans Conjecture 

The most basic arguments surrounding the Jeans conjecture is that 
it is not entirely clear what scale one should average over in order 
to derive the physical conditions in the cloud. In this work, we as- 
sume that the IMF scales with global cloud properties, on scales of 
~ 70 pc. If one were to choose significantly smaller scales, as an 
example, it is possible that the temperature and density would not 
vary in a manner that kept the Jeans mass constant. A number of 
works have posited that a reasonable choice f or scale may b e one 
where the isothermality of the cloud is broken. iLarsorJ 12003 sug- 
gested that this is when the gas becomes energetically coupled to 
the dust, though other works have argued that the isothermality may 
be broken due to radiative feedback (Bate 2009. ; .Krum holz 201 1). 

While we cannot resolve these scales, what likely matters most 
is not the physical conditions in GMCs as a whole, but rather the 
clumps that form stellar clusters. Our model follows the typical en- 
hancements and decrements of the Jeans mass on cloud- wide scales 
given the physical conditions in a galaxy. So long as the physical 
conditions in star-forming clumps within a GMC scale in a self- 
similar manner with the global physical properties of the molecular 
cloud, then the galactic environment may indeed play a role in set- 
ting the fragmentation scale that impacts the distribution of core 
masses. 

From a simulation standpoint, some tentative evidence exists 
that suggests a cl ose connection betwee n th e cloud Jeans masses 
and IMF Work bv lBate & Bonnelll ( l2005h and lKIessen et alj ( l2007h 
showed that the mass spectrum of collapsed objects varies with 
the Jeans properties of the simulated cloud, and that conditions 
comparable to starbursts may have turnover masses as large as 
Mc ~ 15 Mq. This said, the degree of fragmentation is depen- 
dent on the exact form of the equation of state, and the scaling 
of the turnover mass may depend on the numerical resolution for 
isothermal simulations. 

Perhaps the most convincing observational evidence for the 
Jeans conjecture is that it is reasonably well-established that the 
stellar IMF c orresponds well w i th the cloud core mass function in 
some GMCs jAn dre et al." 1996": 'John stone et alJl2doil;lMotte et alj 
l200lHNutter & W ard-Thompson 2007l: lAndre et al.ll201lh . If these 
cores are th e resul t of clump fragmentation that l ead to star forma - 
tion (Evanj ll999l : ITada & Lada 2003; Kennicutt & Evansl l2012l) . 
then a relationship between the characteristic fragmentation scale 
and the IMF is natural. There is additionally some evidence from 
the Taurus star-forming region and Orion Nebula Cluster that sug- 
gest variations in the IMF with clou d density in a maimer that is 
consi s tent with a Jean s scaling (Bricefio et al Jl 20021 ; iLuhman et alj 
I2OO3I ; iLuhmad |2004| ; iDa Rio et alj l2012h . though the dynamic 
range of physical conditions on global cloud scales in the Galaxy is 
limited enough that unambiguous evidence for an IMF that varies 
with cloud Jeans masses is not yet present. ALMA observations 
will allow for better measurements of the relationship between 
GMC properties and core mass functions over a range of physi- 
cal conditions which may shed light on the potential relationship 
between the stellar IMF and cloud Jeans mass. 

Whether the IMF varies with the Jeans properties of clouds, 
or at all with gas physical conditions remains an open question. 
While some observations show mild variatio ns, unambiguous ev - 
idence for systematic variations are lacking dBastian et alj|201oh . 
The Jeans conjecture equally enjoys its share of tentative observa- 
tional evidence, as well as contends against cogent theoretical ar- 
guments. In this work, we remain agnostic as to the validity of the 
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Jeans conjecture, but rather simply aim to understand tlie potential 
consequences of an IMF that scales with global cloud properties. 



5 SUMMARY 

Utilising a combination of hydrodynamic galaxy evolution simu- 
lations, a radiative model for the interstellar medium, and cosmo- 
logical models for the star formation histories of galaxies, we in- 
vestigate the cosmic evolution of the stellar initial mass fimction 
under the assumption that the IMF varies with the Jeans masses of 
molecular clouds in galaxies. 

The IMF can be both top heavy and bottom heavy (at different 
times) for massive galaxies. At early times, the ISM is warm owing 
to the heating of dust by stars, and energy exchange between gas 
and dust. This drives Jeans masses in molecular clouds up, and, 
under the Jeans conjecture, results in a top-heavy IMF. 

As galaxies evolve toward present epoch, their SFRs drop be- 
low Milky Way values, and a lack of cosmic rays allows the temper- 
ature of the ISM to cool below the ~ 10 K characteristic of Galac- 
tic clouds. At these late times, the typical Jeans masses of clouds 
decreases, and the IMF tends toward an excess of low-mass stars 
(bottom-heavy). More massive galaxies quench at earlier times, 
and thus go through more protracted periods of very low-level star 
formation; hence, more massive galaxies at present-epoch have in- 
creasingly bottom-heavy IMFs in our model. 

Because the massive stars that form during the early, top- 
heavy phases have all disappeared by present epoch, the resultant 
M/L ratio for massive ellipticals increases with increasing galaxy 
mass. We find that given plausible galaxy star formation histories, a 
Jeans model for IMF variations provides a reasonable match to ob- 
served M/L — a relations for observed elliptical galaxies at 2: ~ 0. 

In short, by exploring the ansatz that the stellar IMF varies 
with the Jeans properties of clouds galaxies, we conclude that a 
Jeans model for IMF variations in galaxies is able to simultane- 
ously produce bottom-heavy elliptical galaxies at late times, while 
motivating top-heavy IMFs at early times that may help alleviate 
some tensions in high-2 galaxy evolution. We remain agnostic as 
to the underlying Jeans assumption, as well as regarding inferred 
IMF variations at both low and high-z. Improving constraints on 
the IMF in starbursting systems and confirming preliminary trends 
of varying IMF with star cluster age would seem to provide vi- 
able paths towards discriminating the Jeans conjecture model ver- 
sus models that invoke bottom-heavy IMFs throughout the lifetime 
of today's massive galaxies. If the inferred trends suggesting more 
bottom-heavy IMFs with increasingly massive ellipticals alongside 
top-heavy IMFs in heavily star-forming systems indeed are robust, 
then an IMF that varies with the Jeans properties of molecular 
clouds may provide a reasonable physical model to explain such 
trends. 
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